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ABSTRACT 


The Kohn variational method of calculatina scattering phase shifts 
has been examined previously, and calculations of the phase shifts of 
S-wave electrons on atomic hydrogen done, by Schwartz. This paper 
presents the results of using a trial function somewhat simpler than 
Schwartz's and compares the two sets of results. In addition this 
paper presents and compares the results of two different normalizations 
of the variational principle and gives a somewhat more detailed 
development of the variational principle used than in most previous 


papers. 
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I. INTRODUCTION 


In 1961 Schwartz used Kohn's variational principle to calculate 
the scattering phase shift of S-wave (= 0) electrons incident on 


atomic hydrogen. The trial wave function he used had the form 


“Ke (4th, ) aa m fe 


jaws ¢ a San Cg 24" Be) fen /e 


wore $= (1 B,)2 @ EERE + tan) SHER 


G6 PD) fangs 


with as many as 50 terms in the sum [Ref. 1]. 

This paper presents the results of calculations for the same 
problem using a different trial wave function. Assuming the results 
of Schwartz's calculations are correct, a comparison of the new results 
to his is made to determine the relative merits of the new trial wave 
function. The new calculations were also done using two different 
normalizations of the variational principle, and a comparison between 
these results is also presented. As a lead-in to the actual calcula- 
tions this paper first presents a development of the theory behind 


them. 


II. THEORY 


A. VARIATIONAL THEORY 
1. General Form 
First in the presentation of the theory is a general look at 
the variational method of determining the exact wave function of a 


system. The general variational principle starts with the integral 


[= ae | ¥CE-A)Y dT which equals 0 when Y is the 


exact wave function. After some manipulation it is found (See Appendix 


A-1), in spherical coordinates, that 
sl = 242 [SY(E-A)Y dT + lim | PY ESY da 
roo "SS 
— lim | S457 Y dn 


r— ge 
Ci eal Gaull woe Pj or vesy da ~ f ese do 


cr ——} eo 
when V is the exact solution of HY= EY. 





2. Specific Forms 
Now, by specifying the radial portion of Yat r = 0 and as 
roo the general result above yields a particular variational prin- 
ciple. For example, the necessary conditions imposed on Y for a bound 


state of angular momentum £ are that at r= 0 G(r) = OQwhere Y = 





oo b(e,¢) and that as roo 00. These give $I = 0 which 
is the Rayleigh-Ritz variational principle for bound states. The 

Kohn variational method for scattering phase shifts is derived from 
this general formula by putting in a set of conditions appropriate to 


the scattering problem. Therefore an examination of the appropriate 


conditions is now presented. 


B. SCATTERING THEORY 
1. Phase Shifts and Partial Wave Analysis 

First note that the phase shift being calculated is the phase 
difference between the radial portions of the asymptotic forms of the 
wave functions with the scattering potential and without the scattering 
potential. Dividing the incoming wave into partial waves by angular 
momentum (partial wave analysis) there exists a phase shift for each 
partial wave. This set of phase shifts completely determines the 
scattering cross section. In the case of low energy scattering al] 
the phase shifts, except the S or £ = 0 wave phase shift, are very 
small, and therefore the S-wave phase shift is often the only one 
needed for accurate calculations. 

2. Asymptotic Conditions 

So the conditions needed are those on the radial portion of 
the solution for angular momentum & of (E - H)Y = 0 in spherical 
coordinates where H is the Hamiltonian with the scattering potential 
included. At r= 0 the radial portion Ry(r) = Ale) must be finite 
which implies that Xe must equal 0. As re, and the potential 
approaches 0, Y must become like the wave function of a free particle, 
and the most general form X, can assume in spherical! coordinates is 
shown by Schiff [Ref. 2] to be Aylr) rade Elm) Sin(Rr -— +%) 
where 2 is the phase shift of the gt partial wave, and F(%,) is 


an arbitrary amplitude function. 


C. FINAL VARIATIONAL PRINCIPLE 
1. Final Form and Normalization 
Putting these conditions into the general form and allowing 


Ne to be the variable parameter for OT yields the specific form 


Cl. = - $°(%)/k S (See Appendix A-1). 
Now choose F(%,) = secant (7%) and write F (%p) ={ £m) dy 
so that £*(%,) o7 ‘i SF (7) and $I= OF iy, )/k = a Then 
Coal 5 F(m,)/k) = seu + tan(m)/R)=0 which is Kohn's variational 
principle. 

Ordinarily Kohn's principle is written [tan(7)/k) =[+ taal) 
which means that when Y is exact Lt tan(m)/R = tanlm,)/R So by 
taking the variation of L+ meat ge) with respect to several variable 
parameters in a trial function , one of which 1S Ms and finding the 
values of these parameters which make L + tant Stationary, 

§(L + ten (%)/k ) = 0, one gets first order values of 7%» and the 
other parameters. These first order values are then substituted back 


into L+ santa) to give the stationary result Ltan(m)/k], 


2. Usefulness of Variational Theor 





What good is all this? It means that to find the scattering 
phase shifts all one has to do is select a reasonable trial wave 
function with several variable parameters which satisfies the above 
conditions, put it in the integrals, find parameter values that satisfy 
§(I+tan (m/e ) = 0 and put these values back into [tan (ne)/R = 

Lt tan(m)/k If one can find a relatively simple and easy-to-integrate 
wave function that gives good results for a difficult problem, this 


method may save considerable work. This is a big if with complications. 


IJ]. GALCULATIONS 


A. GENERAL OUTLINE 

In the case of variational calculations of scattering phase shifts 
a unique problem appears. In scattering calculations E - H has a 
continuum of eigenvalues which passes through 0. In the calculations 
E - H is represented by a matrix Mi in the subspace of the complete 
Hilbert space of H which 1s spanned by the N terms in the trial function. 
This matrix has N eigenvalues and when any of these eigenvalues becomes 
very small tan(%)/R may become very large and inaccurate. Thus at a 
given energy E convergence to better values with more terms is not 
guaranteed. However, by adding a non-linear variable K to the trial 
wave function and varying it for each value of E the behavior of 
tan(m, )/R can be mapped out. Then the best value of tan(%)/R is 
found by picking the value from the flattest part of the graph of 
tan(%)/R versus K, As the trial wave function becomes more and more 
accurate with more and more variable parameters the flat regions of 
the curve become flatter and longer; converging, hopefully, to better 
and better values of the phase shift. When the improvement becomes 
small compared to the size of tan(%,)/k the final value is determined 
from the flat region of the graph and the probable further improvement 
due to more terms. The error is estimated from reading the graph and 
how well the further improvement can be determined. A good presenta- 
tion of this general method of calculation is given in Schwartz [Ref. 3]. 
‘ The method used in the present case differs by using the values from 
the curves directly without estimating probable improvement and 


estimating the error only from reading the graph. 
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B. THE HAMILTONIAN 

Now apply the above to scattering S-wave electrons off atomic 
hydrogen. In this problem there is an incoming electron, the hydrogen 
proton and the electron bound to the proton, all of which interact with 
one another. To take in all these interactions the following 
Hamiltonian was used 

22 E-H]= ok - +r t+ - + Ma Fe 
where radial distance and energy are in Bohr radii and Rydberg units 
respectively (See Appendix B). Therefore Ro is the energy of the 
incoming electron, and -1 is the energy of the bound electron in the 
ground state. To avoid the complications of excitation of the bound 


electron R was limited to less than .866, which corresponds to about 


10.2 ev, so only low energy scattering was considered. 


C. THE TRIAL WAVE FUNCTION 
However one still had to take into account the indistinguishabi lity 
of the electrons, and this was done with the trial wave function. The 


following wave function was used 


-K(yitn), . i 
Cjme | (ee te) are 


(\ 
Ps 
+ 
4 
na 


Ss 


> = ge (12P,)[ jolkr) - tanta) j, (ke) _"* 


where Be is the parity operator that interchanges ry and ro and where 





Sin RG Sinky, cosRr. 
Jo (Rr) = Ry, and Jy (Rr, ) > Teel = are spher- 


ical Bessel functions. The + sign gives a symmetric wave function which 


implies spin angular momentum S = 0, and the - sign gives an anti- 
symmetric wave function which implies S = 1. Thus the 1 - Rio takes 
care of the indistinguishability of the electrons and their fermion 
properties. 


is the asymptotic part of the wave function and approaches 


ve Sin(krt+%) .-%G 
air secant (7) _. © as r)—? © 


\ 


and “mT 


JE secant (m) sin eeegl e™ 
Ps 


S oo 
a ro 


contains the necessary asymptotic plane wave condition for one 
electron, coupled with a ground state condition for the other electron 
in the e ’ term. The sum portion contains the close-in terms with the 


variable parameters CoH Note that # is also a variable parameter. 


D. INTEGRALS AND FINAL FORMULAS 

Now with two electrons two radial coordinates, ry and ro» and 
their difference, Tyo» are used. Does this change the variational 
principle? The answer is no, with the proper normalization of the 


trial wave function (See Appendix A-2). So the variational principle 


remains [tan (ne) /R] = [+ tant, )/R . Now defining 


G+ ECi Xin, A= tan(y) ) C.F Oe BIG: 





22 (x (e-Ald drdt.= Rit aR: 


Miss S@ | Xi (E-A)X) eT dT 
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one solves lila A/R ) = ra Mier Goat aae and 


raf 2 
JA et Aff) =O to get the C. and A to first order. Then 


inserting these values back into the variational equation gives 


2 
p w,+B, t+ I/k 
[a/R] = mies ei 


gaerreiera() .() Sa) . See 
where Wy = = C-R., W, 27-C5R: and W, 2 C:R, as shown by 
Schwartz [Ref. 3] . These are the final integrals and equations used 
in the calculation of the phase shift using F(%) 2 secant(y). The 


phase shift was also calculated using f(7) > cosecant (7) which gives 


[- cot in)/R ] = J - cot (») /k and requires 


¢= HE (12 P,) Loot) jeri) ~ jy lawl} e > 
The final formula is 


an 
: w, + B,> I/k 
[- a/R] = Wat Ba - BAL 


where A= cot (x) and the W's and B's remain the same as before except 


Lae 5 0 ] en) ] 
B= A Bayt AaB Bee AC. #C., and R, AR. + Ra. 


E. THE ACTUAL CALCULATIONS 
The actual calculations are done as follows. First the integrals 


By: Ban « s IRE are done by hand for the general case, and then 


the final numbers are obtained by computer evaluation of the general 


integrals for particular cases. Then 2M 5¢ jamie is solved independ- 
ently of the first order solution for ¥ since a co FAC, and Mi 

(9 = ] ' 
independent of A implies ya es rams RY and EN; 5. j = -R.. es 
solution is most accurately done by matrix inversion. The (Z| 


factors in the close-in terms of Y are there to increase the accuracy 
of the matrix inversion and have no effect on the value of the phase 
shift. Finally the W's and the stationary value of 2 are calculated 
using the previously given formulas. 

The above calculations were done for values of K from 0.6 to 2.5 
in steps of 0.1 for each value of R from 0.1 to 0.8 in Steps of 0.1] 
with up to 40 close-in terms. This was done for both the symmetric 
and antisymmetric wave functions and both secant and cosecant normal- 
izations. Then the values of tan(%)/R or -cot(7) /k were plotted 
versus kK for each case by the computer (See Figs. 1, 2, 3, 4), and 
values of tan()/k and -cot(n) A selected from the flattest portion 
of each curve. The actual values of 7 were then calculated from these 
numbers for comparison with Schwartz's values. Additional values were 
obtained for R= 0.4 at 49, 56, and 72 terms to see how much better 
the values might get and where significant inaccuracies in the matrix 


inversion would occur. 
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Fig. 1 - Graph of tan(»)/k versus kK (kappa) for R=0.5 and S = 1. 
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Fig. 2 - Graph of tan(7)/R versus K (kappa) for R= 0.5 and'S = 0. 
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Fig. 3 - Graph of -cot(m)/k 
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versus Kk (kappa) for R= 0.5 and S = 1. 
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Fig. 4 - Graph of -cot(y)/k versus K (kappa) for R= 0.5 and S=0. 
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IV.) DIStuSsion OF RESUEIS 


A. COMPARISON WITH SCHWARTZ'S RESULTS 

Table I lists the values of tan(#)/R and 7 obtained from the 
present calculations and from Schwartz's calculations for S = 1. Table 
II does the same for S = 0 and Tables III and IV do the same for 
-cot(7) /k for S = 1 and S = 0 respectively. The uncertainties in the 
present values listed in the tables represent the errors one would 
estimate from the graphs in the absence of Schwartz's values for compar- 
ison. Comparison of the results 1s done quite easily from the tables, 
and the following results are obtained. In Table I the percentage 
differences between the values of present calculations for 7 and 
Schwartz's calculations are between 2% and 8% with one case of 12%. 
In Table II the differences in are between 10% and 20% with one 
case of 40% and one of 207%. In Table III the differences in 7? are 
between 2% and 8% with one case of 15%, and in Table IV the differ- 
ences are between 9% and 20% with one case of 30% and one of 213%. In 
the cases with 207% and 213% differences in” the differences in the 
absolute values are only 2% and 9% respectively but the values have 
opposite signs. These two cases are the only such cases of difference 
in sian, both occur at R= 0.3 for S = 0, and neither can be explained 
adequately at the present time. Ignoring these two very large cases, 
the average differences from Schwartz's values in ® for both secant 


and cosecant normalizations, are 5.5% for S = 1 and 15.9% for S = QO. 


We 


Table I. List of present secant normalization values versus Schwartz's 
values for S = 1. The numbers in parentheses give the uncertainty in 
the last digit. 


se 
ti 


-2.460(5) -0. -0.2028(4) 





















=28555i( -0.4245(5) 


-0.6419(8) 










-3.125(5 


5) 
-2.745(5) 
) -0.84776(35) 

) 


-3.723(8 -1.0370(4) 






-5.85(6) -1.2087(7) 










-8,35(2) -1.3619(6) 





-1.498(3) 





-28.75(5) 


Table II. List of present secant normalization values versus Schwartz's 


values for S = 0. The numbers in parentheses give the uncertainty in 
the last digit. 


fn ome | —_lreion) | tote | pasos 


0.1] -10.00(5) -0.785(2) -6.68(2) -0.589(1) 
-1.276(5) -9.23(2) .0743(9) 
e465) -26.4(1) .4452(5) 
ee 2) 15.87(4) .4145(4) 
1.062(1) Sel 2) .2017(13) 
0.908(1) 2.845(8) .041(1) 
0.7946(3) ¥.917(5) .9303(12) 
Ge72 193) y.530(5) -8857(16) 





Table III. List of present cosecant normalization values versus 
Schwartz's values for S = 1. The numbers in parentheses give the 
uncertainty in the last digit. 


-0.23417(3) -0.2028(4) 



















.4597(12) .4245(5) 





a58023(.1) 






.6419(8) 





.883859(5) 





.84776(35) 





1.09847(3) .068541(11) .0370(4) 









.2399641(27) .2087(7) 






0.572425(5) 






07125708(6) .39275(4) .3619(6) 








0.056993(7) 52523316) .498(3) 





Table IV. List of present cosecant normalization values versus 


Schwartz's values for S = 0. The numbers in parentheses give the 
uncertainty in the last digit. 


ih 
Ol : : 


ize sire -0.72686(13) 
1.6855( 2) -1.24566(4) 
-0.1188(1) Was3517(3) 
-0.77845(2) 1.268932(7) 









= 
NO 





















-1.1076(2) .06504(7) 








| -1.29527(3) 2 10137(11) 













-1.39933(3) .79574(1) .9303(12) 












-1.409835(5) .725378(2) .8857(16) 
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B. COMPARISON OF DIFFERENT NORMALI ZATIONS 

The average difference between the values of 7 calculated with 
secant and cosecant normalizations is 2%. One other observation is 
worthy of note, and that is that the values of -cot(%)/k found by the 
cosecant normalization have more apparent significant figures than the 
values of tan(y)/k . This is due to the fact that the plots of -cot(%) /k 
versus ‘ct do not show fluctuations in -cot(7)/k anywhere as large as 
in tania) k . To see this compare Fig. 1] with Fig. 3. Since the 
normalization of the variational principle is arbitrary there should 
be no difference in accuracy from one to the other, and the comparison 
with Schwartz's results shows this to be true. Without Schwartz's 
results for comparison one would seem to have more confidence in the 
results of the cosecant normalization due to the larger number of 
apparent significant figures. However, without a standard of compari- 
son, there is no reason why one normalization could not give more 
Significant figures in a given case and still not give any better an 


answer for 


C. ACCURACY AND CONVERGENCE 

The results examined above are for 40 terms in the sum and they show 
fairly good results for S = 1 and generally poorer results for S = 0. 
An examination of runs at R=0.4 with 49, 56, and 72 terms shows the 
following. Beyond 40 convergence of the flat region toward Schwartz's 
values is very slow and not always guaranteed, since for the S = 0 
‘secant normalization at 49 the value actually became worse. In addition, 


at 56 and 72 terms, significant errors in the matrix inversion were 
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noticed ana the accuracy beyond about 50 terms became questionable. 
The limit imposed by computer accuracy was therefore between 50 and 
60 terms, and within this limit continued convergence did not give 


appreciably better results. 


D. RELATIVE ACCURACY OF S = 1 AND S= 0 

The triplet (S = 1) calculation is better behaved than the singlet 
(S = 0) calculation, probably for the following reason: the anti- 
symmetric (triplet) Y is automatically 0 for r) = Py $0 that the 
electrons cannot coalesce. This effect is qualitatively the same as 
the effect of electrostatic repulsion, and thereby operates to remove 
some of the burden from the variational machinery. The singlet wave 
function contains no such built-in convenience. Thus for S = 1 the 
calculations have a greater absolute accuracy than for S = 0 and so 


should have a greater relative accuracy for S = 1] than for S= 0. 
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Vv. CONELUSIONS 


There are two differences between the trial functions used in this 
paper and those used by Schwartz. In the asymptotic part Schwartz 


cos(kr,) 


-Ky/V, 

sed Thee blewe =y instead of Jy (Rr) and in the close- 

in terms Schwartz included an rh which was left out in the present 
calculations to make the integrals easier. From the comparison of the 
present results with Schwartz's results one concludes that in general 
the convergence of the curves for the present trial function is not 
as good in an equivalent number of terms. The idea that the effect 
on convergence of eliminating the a factor could be offset by using 
more terms of a simpler form seems to be cut short by the inaccuracies 
in the matrix inversion and slower rates of convergence with larger 
numbers of terms. Therefore one concludes that a more complex trial 
function must be used to get good results, and that one can be misled 
by the apparently qood accuracy (the small uncertainties in the tables) 
of the simpler trial functions without comparison. Finally one can 
conclude from the above comparison that the simpler form works 
reasonably well for S = 1 but not so well for S = 0, and that there 

is really no difference in accuracy between the two normalizations 


despite the larger number of apparent significant figures for the 


cosecant normalization. 
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APPENDIX A 


THE VARIATIONAL PRINCIPLE 
This appendix shows the steps in the derivation of the variational 
principle used, for the case of one radial coordinate r, and for the 
case of two radial coordinates ry and ro. 
1. One Radial Coordinate 
ape ( Y(E-R)Y d= where E is a constant 


bh 2 
energy term and H= -tm V+ U(r) 


Vi 


Start with I 


where U(r) is a central potential. Then write 
ST= 42 [Ssvle-A)yar + SY(E-W)s¥ dT 
Now using «= WWSY = D(H VSY) ~ Ve (S4¥0Y) + SHV'Y 


One gets Rie a 2m f5¥(E-A)Y dr + fF (YTSY) dT 


— (FS49W) dr. 


men srs 292 S sate-Alyar + So wh se dS — L544 g8 


Se Psy(e-AluaT tlim [Sor 4S da-Sr or eda] 
ie Aa) 
Since the integrals are over all space and ig arrdin. « 


Finally 


~ id PY, 
ST= lim Ser YSS¥da- lim foo S4 se ¥dn 
Vr 00 I~) 08 
for 'Y exact. With ¥Y properly normalized the integral over the angular 


part can be made unity. Then 


Se him [Re rae SRe ~ SRe PSF Re] 


where Rg is Ae. radial portion of '% Then inserting the condition 
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on Rg at & 


LT 
sin (Rr - 3 +H) 


int =a §(7.) Rr 
one finds, using SRe ? in, Ry 5% 
that § I= -F(ne)/R 5% for exact. Note that in 


the case of one radial coordinate the correct normalization of Y TOT 


L~ 0 is wat 
2. Two Radial Coordinates 


> 


In this case H = ae (Vv, +Vi) T Neel, pe) 


Then as before 


T= lim LSP YES dadte - limdule 54S Yda, dh 


Y.- 0c A 
tml n Ys SY dadh - Limb S1ka Yalan; 
+ 0 


for ¥Y exact. Aaain with Y properly eewaT i zed the integrals over dnd. 
and dQ.dQ, can be made unity so 
el} ve am OF 
sI= lin SL Rese SRe- r*SRede Rede Sr tlinl, LA wierd SRe- Fd Kear Ryliedr, 
G0 


Then inserting the conditions 


Esc ta) 8 = in(ke- Shay) 
Bere 60 “e S ; Rg aon F(y,) sina tel 


one gets §T= of ul oy. where the 4—can be eliminated 
by putting an additional factor of = in the normalization of ¥. So 
finatiy S$ L= 7 £(nel/k sNe as before. Note that the 


: ; : Pa 
correct normalization of for L= 0 is now ——. 


cart =. 


APPENDIX B 


THE HAMILTONIAN 


t rs ae 3 ' 
Write GY = Sg tT Yt e7( a ete) Y= EY 
(2 ue 
ge a | " tk ag Va Se 
then Toe +Vi eee - ara + Y= ( am 2 2) 
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